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We define the flatness and quasi-flatness problems in cosmological models. We seek solutions to 
both problems in homogeneous and isotropic Brans-Dicke cosmologies with varying speed of light. 
We formulate this theory and find perturbative, non-perturbative, and asymptotic solutions using 
both numerical and analytical methods. For a particular range of variations of the speed of light 
the fiatness problem can be solved. Under other conditions there exists a late-time attractor with 
a constant value of Q that is smaller than, but of order, unity. Thus these theories may solve the 
Q.^ ] quasi-fiatness problem, a considerably more challenging problem than the fiatness problem. We also 

. discuss the related A and quasi- A problem in these theories. We conclude with an appraisal of the 

' difficulties these theories may face. 

§ ■ I. THE FLATNESS AND THE QUASI-FLATNESS PROBLEMS 

^ ! 

i The observable universe is close to a flat Friedmann model, the so-called Einstein-de Sitter universe, in which the 
energy density, p, takes the critical value, pc, and the homogeneous spatial surfaces are Euclidean. All astronomical 
evidence shows that we are quite close to this state of flatness, although a value of flo in the vicinity of 0.2 is preferred 
by several observations. 

' It is therefore disquieting to notice that the flat Friedmann model containing dust or blackbody radiation is unstable 
, as time increases. Small deviations from the exact $1 = 1 model grow quickly in time, typically like a^, where a is the 
T-H ■ expansion factor of the universe. The observed fio ~ 1 state therefore requires extreme fine tuning of the cosmological 
, initial conditions, assumed to be set at Planck epoch, because a has increased by a factor of order 10'^^ from that 
epoch to the present. This is the flatness problem. If the universe is slightly open at Planck time, within a few Planck 
times it would become totally curvature dominated. If it is initially slightly closed, it would quickly collapse to Planck 
f~i ■ density again. Explaining its current state requires an extraordinarily close proximity to perfect flatness initially, 
or some sequence of events which can subsequently reverse expectations and render the flat solution asymptotically 
Q , stable. 

$^ ' Particle physics theories naturally contain self-interacting scalar matter fields which violate the strong energy 
c/j ,, condition (so the density and pressure, p, obey p + 3p < 0). These can make the flat solution asymptotically stable 
Ci ■ with increasing time, allowing the asymptotic state to naturally be close to flatness. Cosmological histories in which 
J> I a brief period of expansion is dominated by a matter field or other effective stress which violate the strong energy 
' • condition, and so exhibit gravitational repulsion, are called "inflationary" . 

] Solutions to the flatness problem have been proposed in the context of inflationary scenarios ||l|, pre-Big-Bang 
^ models and varying speed of light cosmologies ||-^. In all of these theories, Q — 1 becomes an asymptotically 
stable attractor for the expanding universe. The observed r2o ~ 1 state results then from a temporary period of 
calculable physical processes, rather than from highly tuned initial conditions. In such scenarios the price to be paid 
is that flo should be very close to unity, |ilo — 1| < 10~^, if one is not to invoke an unmotivated fine tuning of the 
initial conditions again. 

If we take the trend of the observational data seriously, then explaining a current value of Qq of, say, 0.2 — 0.9 is 
yet another challenge. We call it the quasi-flatness problem. Solutions to this problem have been proposed in the 



> 



context of open inflationary models |lj 
Anthropic Principle pq] is usually invo 



In these one has to come to grips with some degree of fine tuning. The 
ced for this purpose but considerable uncertainties exist in the range of 
predictions that emerge and there does not appear to be scope for a very precise prediction of f2o in, say, the range 
0.2. Undoubtedly, it would be better if one could find mechanisms which would produce a definite of order one, 
but different from 1, as an attractor. In a recent letter we displayed one theory in which this possible. Here we 
present further solutions in support of such models. We explore analytical and numerical solutions to Brans-Dicke 
(BD) cosmologies with a varying speed of light (VSL). These generalise our earlier investigations of this theory 
We show that if the speed of light evolves as c (x a" with < ti < — 1 there is a late-time attractor a.t — —n. Hence, 
these cosmologies can solve the quasi-flatness problem. This work expands considerably the set of solutions presented 
in ref. B- 
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We note that the existence of the dimensionless fine structure constants allows these varying-c theories to be 
transformed, by a change of units, into theories with constant c, but with a varying electron change, e, or dielectric 
'constant' of the vacuum. This process is described in detail in ref. 0] where a particular theory is derived from an 
action principle. A different varying-e theory has been formulated by Bekenstein [pl but it is explicitly constructed 
to produce no changes in cosmological evolution. A study of this theory will be given elsewhere ]T^ . 

There also exist analogues of the flatness and quasi-flatness problems with regard to the cosmological constant. 
The lambda problem is to understand why the cosmological constant term, Ac^ in the Friedmann equation does not 
overwhelmingly dominate the density and curvature terms today (quantum gravity theories suggest that it ought to 
be about 10^^° times bigger than observations permit |^,^). The quasi-lambda problem is to understand how it 
could be that this contribution to the Friedmann equation could be non-zero and of the same order as the density or 
curvature terms (as some recent supernovae observations suggest pO|). 

In Section || we write down the evolution equations in these theories. The varying G aspect of the theory will 
be accommodated in the standard way by means of a Brans-Dicke theory of gravitation. This theory is adapted 
in a well-defined way to incorporate varying c. We write equations in both the Jordan and Einstein's frames. In 
Section III we study solutions to the flatness problem in the perturbative regime (SI « 1), when both c and G may 
change. In Section IV we present non-perturbative solutions when G is constant, and in Section ^ when G is varying. 
In Section VI we derive simple conditions on any power-law variation of c with scale fac tor if the fla tness problem 
is to be solved. These conclusions are reinforced by some exact solutions in Section VII. In section VIII we show 
that the quasi-flatness problem can be naturally solved in a class of varying-c cosmologies and then we discuss the 
solution of the quasi-lambda problem by such cosmologies in Section IX. We conclude with a summary of our results, 
highlighting the cases in which we can claim to have solved the quasi-flatness and quasi-lambda problems. 



II. COSMOLOGICAL FIELD EQUATIONS 

In BD varying speed of light (VSL) theories the Friedmann equations are 1^ : 
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a\ _8Trp Kc^ 
a J 3(l> 




(2) 



where the speed of light, c, is now an arbitrary function of time, K is the curvature constant, and a; is the constant 
BD parameter. The wave equation for the BD scalar field (p— l/G is 

^-)-3-<^=-^(p-3p/c2) (3) 
a 3 + 2lj 

From these three equations we can obtain the generalised conservation equation. Since after we impose the equation 
of state, the varying c term only appears in eq. (H), the only new contribution to 2da is from the —2Kcc term. Hence, 
together, these imply the 'non-conservation' equation: 

In the radiation-dominated epoch {p = pc^ /i) the general solution for </> is 

% (5) 

If the integration constant a = the usual solutions for VSL in general relativity with varying c follow illustrating the 
fact that for a radiation source any solution of general relativity is a particular solution of BD theory with constant (j). 
In the next sections we explore the solutions to the fiatness problem when (/i 7^ 0. We will also consider the radiation 
to matter transition in these theories. 

Equations (|l|)-(^ a-pply in the so-called Jordan frame. It will be useful to introduce the Einstein frame, by means 
of the transformations. 
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di = y^dt (6) 

d = (7) 

(T= (tJ + 3/2)^/^ln(G0) (8) 

p=(G</.)-V (9) 

p = (G(/>)-2p, (10) 



which are performed at constant c. These may be regarded as merely mathematical transformations of variables. 
The Friedmann equations in the new frame are 

5—^^^[p + 3p/c']-— (11) 
a 3 6 



d'Y Kc^ SnGp 



d J d^ 3 6 

where ' = d/dt. The transformed scalar field equation for a is 



(12) 



a" + 3^a' = -^={p-3p/c') (13) 

These are just the standard Friedmann equations with constant G and a scalar field added to the normal matter. The 
scalar field behaves like a 'stiff' perfect fluid with equation of state 

In the Einstein frame, if c = 0, the total stress-energy tensor is conserved, but the scalar field and normal matter 
exchange energy according to: 

p' + 3^{p+p/c')^-{p'^ + 3^{p.+p./c')) = ^f^'^ (15) 

a a V6 + 4a; 

If c 7^ 0, one has instead 

p' + i-r{p + P c^)^ ^^_^ -^ + -— ^- 16 

V , oO'/- , - / 2n Cr'{p-Sp/C^) 

Pa+^—{Po+Pa C)= (17) 

a V6 + 4w 

All equations derived for standard VSL theory, with constant G, are valid in the Einstein frame. However, one should 
always remember to add to normal matter the scalar field energy and pressure (so the total density and pressure are 
given hy pt ^ p + Per and pt ^ p + Pa, respectively). 



III. PERTURBATIVE SOLUTIONS TO THE FLATNESS PROBLEM 

We first study solutions to the flatness problem when there are small deviations from flatness. Let us define the 
critical density, , in a B-D universe by means of the equation: 



d \ Snpc 4>a u cf) \ 



J 30 0a 6 I 
that is, Eqn. with K = 0. In the Einstein frame the critical density in normal matter is 



(18) 
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In terms of total energy density pt = p + Pa , the critical energy density in the Einstein frame is: 

3 /a 



Accordingly, we may define a relative flatness parameter: 

et = (21) 

Pet 

As shown in the previous section, the usual equations for VSL theory should apply to this quantity. Therefore, one 
has 

= (l + e()et(37-2)^ + 2-et (22) 
a c 

with the equation of state given by 

7 - 1 = p/pc^ - p/pc^ (23) 

with 7 constant. 
But, since 

et = T^, (24) 

i + (2'- + 3)s4; 

we see that the natural quantity to quantify deviations from flatness in the Jordan frame is not 

Pc 



but an adaptation of it, 



which satisfies the equation 



^ a 2 

If 5 ^ 1 this can be integrated to give 



(25) 



S = P < e, (26) 



(5= (l + ,5),5(37-2)( - + ) +2-5. (27) 



5 cx a^-'-^c^^^ (28) 



Hence, we have the solution 

1 C (2w + 3)02 



(29) 



If G oc varies then the second term on the right always works against solving the flatness problem. Hence solving 
the flatness problem in BD requires that </> decreases in the early universe and that the first term on the right of eq. 
( p9| ) dominates the second. 
In a radiation-dominated universe (7 = 4/3): 



1 _ C (2w + 3)a2 



(30) 



where a was defined by eqn.(|^). If a is small and c = cqo" one still solves the flatness problem when n < — 1. Notice 
that since pc oc Xja^ at late times, the second term always eventually becomes negligible. 
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FIG. 1. The evolution of e from various initial values for cosmologies with (p = const and c = cq. We see that the flat universe 
is unstable. 



IV. NON-PERTURBATIVE SOLUTIONS WITH = CONST 

A. Matter and radiation-dominated cases 

In order to explore non-perturbative solutions to the flatness problem we solved Eqns. (0)-(^ numerically. First 
we consider the case where (/> = (j)Q. Exact solutions exist in this case for matter and radiation-dominated universes. 
Setting (j) — 1/G in Eqn. ( p7|) leads to 

e= (l + e)e(37-2)- + 2-e (31) 
a c 

which, assuming c/c — na/a, becomes 

e = e((l + e)(37-2) + 2n)-. (32) 

a 

The general structure of the attractors can be inferred from this equation. However, it can also be integrated to give 



1 



2n 
37-2 



1 - Ca-(2"+37-2) ' 



(33) 



with n < 0, and where the integration constant C is chosen so as to enforce the specified initial conditions. Our 
numerical code matched this analytical result to very high accuracy. In Fig. |^ we plot numerical results for a 
radiation-dominated universe. We comment on four possible situations, illustrated in Figs. |]-^. 

If c is a constant (n = 0) then the flat universe (e = 0) is unstable (Fig. |^). There are two attractors at e = — 1 
and e = oo, showing that the universe tends to become curvature dominated. For slightly closed universes this means 
evolution towards a Big Crunch singularity. For slightly open universes it means evolution towards a Milne universe 
(e ~ —1) which in the case of constant G is simply an expanding empty spacetime with a oc t. If the speed of light 
were to increase with a (c cx a" with n > 0) the situation is the same, but with a stronger repulsion from e = 0. 

If c oc a" with n < —(87 — 2)/2, then, as pointed out in Jp|, the flatness problem is solved if e is not too far from 
zero initially. In fact, e = is now an attractor (see Fig. 0). The non-perturbative analysis reveals two novelties: 
there is an unstable node at e = —1 — 2n/{Sj — 2) and e = 00 is an attractor. This means that if the universe at any 
given initial time satisfies e > — 1 — 2n/(37 — 2) it will not evolve towards the flat attractor. Instead, it will evolve 
towards a big crunch with e = 00. Therefore, for closed universes, the fiatness problem is solved only if they are not 
too positively curved. 

It is curious to note that the situation is different for open universes. No matter how far they are from fiatness they 
will always evolve towards the Einstein-de Sitter model at late times. This was first pointed out in j^, where it was 
noted that in VSL theories the Einstein-de Sitter universe is an attractor even if the universe is initially Minkowski 
(empty) or de Sitter (dominated by a cosmological constant). 

The greatest novelty, however, appears for the case in which —(37— 2)/2 < n < 0. Although e = is now unstable, 
there is a stable attractor at e = — 1 — 2n/(37 — 2) (see Fig. ||). The system is now repelled from e = — 1 but is 
attracted to e = 00. 
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FIG. 2. The VSL solution to the flatness problem. We have assumed a radiation-dominated universe, and considered a 
situation with n < —(87 — 2)/2. We see that although the flat universe is now an attractor, if the universe is ever too closed it 
will end in a Big Crunch, rather than approach the flat attractor. 
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FIG. 3. The VSL solution to the quasi-flatness problem. We have assumed a radiation-dominated universe, and considered 
a situation with —(87 — 2)/2 < n < 0. The attractor is now at a value of f2 of order 1, but not necessarily 1. Closed universes 
end up in a Big Crunch, but any open universes are pushed towards Q, = — 2n/(87 — 2). 
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FIG. 4. The transition case n = —(87 — 2)/2 illustrated for a radiation dominated Universe. 
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Hence we have the following scenario: all closed universes evolve towards a big crunch, no matter how close to 
flatness they are initially. Thus a selection process excludes even slightly closed universes. In contrast, flat and all 
open models (no matter how low their density) evolve towards an open model with fl = ~-2n/{3j — 2), a value that 
lies between and 1, and is typically of order unity. We therefore have a solution to the quasi-flatness problem, for 
open, but not for closed universes. 

The n — —(87 — 2)/2 case is a transition case in between the last two situations just described (Fig. U). In this 
case e = is a saddle: stable if approached from below but unstable from above. 



B. The radiation-to-matter transition 



In quasi-flat scenarios the major contribution to the matter density of the universe is the term resulting from 
violations of energy conservation. This can be seen from 

B 3Kcn(j)n 

^-^ + 4.(37 + 2n-2) -- ^''^ 

derived above. Taking radiation as an example (7 = 4/3), if n < — 1 we are pushed towards the attractor fl — 1 
and the energy-violation term (second term) becomes negligible. In scenarios in which the quasi-flatness problem is 
solved the converse happens: the second term dominates. Matter is continuously being created and this provides the 
dominant contribution to cosmic matter at any given time. Matter which was created in the past gets diluted away by 
expansion very quickly. In effect, we have a state of permanent "reheating" . This scenario has a passing resemblance 
with the steady state universe, in which the so called C field continuously creates more and more matter. 

If we couple a changing c say to the Lagrangian of the standard particle physics model then all particles are 
produced, and whether they behave like matter or radiation depends on whether T ^ m or T ^ m. If7i<0 the 
energy density in created matter decreases, and the universe cools down as it expands. Hence there will necessarily 
be a radiation to matter transition in these scenarios, just like in the standard Big Bang theory with constant G and 
c. 

A complication arises when we attempt to model the evolution through this transition because the structure of the 
attractors changes. For a radiation-dominated universe one requires — 1 < n < for the open attractor Q — ~n to be 
achieved. For a matter-dominated universe the attractors for open universes are given hy Q = —2n and are achieved 
when —1/2 < n < 0. Hence, if —1 < n < —1/2 we have attraction towards a (quasi-flat) open universe during 
the radiation-dominated epoch which is then driven towards flatness in the subsequent matter-dominated epoch. If 
— 1/2 < n < we have attraction to a (quasi-flat) open model in both epochs but the universe will evolve from a 
trajectory asymptoting to = —n toward one asymptoting to = —2n after the matter-radiation equality time. The 
solution (RSh is then no longer valid because 7 is time dependent. 



V. SOLUTIONS WITH VARIABLE (/> 



A. Exact Radiation solutions 



For radiation, 7 = 4/3, we can solve the field equations by generalising the method introduced for scalar-tensor 
theories introduced in ref. [|l^. Define the conformal time 

adr — dt (35) 

so ' = d/dr. Eq. (||) integrates to give 

A 



with A constant . 

Now change the variables of eq. (0) by introducing 



and d/dr. It reduces to 



,2 



y = (f>a^ 



(36) 
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2u} 327rj/pa* 



2„,2 



- 4A'c 2/ 



For radiation, eq. (&§) gives 



pa^ = B 



SKcQiny 



47r(2m + 1) 



If we substitute (pi) in then we have 



2lu, 32t:B 



-y + AK4 



2m 



(2m- 



2(m+l) 



(37) 



(38) 



(39) 



This can be integrated exactly for appropriate choices of m and its quahtative behaviour is easy to understand. 
Note that if we can solve ( ^9| ) for y{T) — (f>a^, then we know <p{T,a) and can solve (|3^) to get 0(t) hence a(T), and 
finally a{t) from ( ^ since 

5- = - 40 

We note that the curvature term in ( |39| ) changes sign for a special value of m, 

1 

m = = 

2 

Examining ( |39| ) we see that the = solution appears to be approached if 2{m + 1) < 1, ie m < —1/2. 
When m = — 1 there is a simple special case which allows an exact solution. Eq. (p9) is now 



y'^^a + Py + T 



where a,f3>0,T^ +4,Kcl. We integrate (|oh to get 



P / -.2 a + T 2 



Hence, 



and 



/3(t + ro) - 2^/^TV 



_/3(T + ro) + 2VcHT' 

y ^ P^jr + T^f \P{t + To) + 2x/^TT^ -"^^^ 
4/3(/.o 



/3(t + To) - 2^/^+f 

We see that if a + F > then, as t ^ oo, we have a — > t ^ t^/^ unless the denominator blows up. 
If a + F < then we get 



2A ^,^^-.M(r + ro) 



V-(" + r) V V-(" + r) 



and 



1 ro2/ N x2n r -i//3(t + to) 

4^ f''^" + + exp{-- tan (^^^ 



} 



where 

= -4(a + F) 

Looking back at the defining equation (|39|), we see that for this m — —1 case, F = 4iirc§ and so 

a + F = A^{1 + '^) + AKcl 
and for large A we have a + F > when K > Q. 



(41) 



(42) 
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FIG. 5. The evolution of e = Q — 1 and of = 1/G from a Milne starting point, with 7 = 4/3, c = co/a^, and different values 
of the integration constant a, of eq. (^). We see that flatness is achieved faster with a changing G (increasing or decreasing). 

B. Numerical solutions 

A numerical evolution of the equations with a varying (f> reveals that for reasonable initial values of the structure 
of attractors is not affected. However, the speed at which attractors are reached may be increased if </) is allowed to 
change. As an example we shall consider the case in which one starts from a Milne universe (e = — 1). In |^ it was 
argued that this is the natural initial condition to consider in the context of VSL cosmologies. 

Let us then consider the case 7 = 4/3, and integrate Eqns. (|l| ), (||), and (|^). We consider solutions of the form (||) 
with various integration constants a. In Figures |^ and |^ we show the evolution of e — fl — 1 and cf) from such a state, 
with c = cqu" and n — —2, —1/2. In the first case we see that we still have a flat attractor, which is reached much 
more rapidly if (f> is increasing or decreasing. In the latter case we have a quasi-flat attractor, with e = —1/2, whether 
or not (j) is allowed to change. The attractor is achieved faster with a changing 0, especially a decreasing one. 



VI. ASYMPTOTIC SOLUTIONS TO THE FLATNESS PROBLEM 



Exact solutions are only possible for particular equations of state and c- variation laws but it is possible to understand 
the asymptotic behaviour in general. The VSL theory Brans-Dicke solutions are given for eqns. (0)-(|) together with 
(^. We assume a perfect- fluid equation of state given by (p3|). In the flat case {K = 0) the equations reduce to those 
of standard BD flat universes with constant c. At late times, when K = 0, the general BD solutions for all 7 approach 
the particular ('matter-dominated' or 'Machian') power-law solutions ||l7| 



a{t) cx 



(43) 
(44) 



where 



2 + 2lu(2 - 7) 
4 + 3a;7(2 - 7) 



(45) 



D = 



2(4-37) 



4 + 3u;7(2 - 7) 



(46) 



Hence, we see that D + 3jA = 2; note that = for radiation;. These exact power-law solutions reduce the special 
(/) — const, a oc t^^^ general relativity solution in the radiation case. We shall therefore look at the stability of these 
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FIG. 6. The quasi-flat scenario with VSL and a changing G. We consider a Milne starting point, with 7 = 4/3, c — co/a^^"^ , 
and different integration constants a. The quasi-flat attractor is still the same as in the constant case, but it is achieved 
faster with a changing (f> (especially a decreasing one). 



K = Q asymptotes when we turn on the K ^ terms in eqns. (|)-(|) and d) with c{t) included. If we substitute 
(Es) and (Hq) in (0) with c = cqo" then we have 



+2A(n-l) 



47r 



47r 



Integrating, we get 



pt 



37A 



B 



87r(l + An~ A) 



^2{l+An~A) 



(47) 



if A{n — 1) ^ —1. Thus, to solve the flatness problem we will need the B term to dominate the K term on the 
right-hand side of eq. (|47|) at large t; that is, since A > for expanding universes, asymptotic approach to flatness 
requires 



n < 1 



1 

A' 



Using (Esh, this gives the condition 



n < 



(2-7)(2-37)c^-2 



2[l + c^(2-7)] 

For radiation and dust universes this condition for solving the flatness problem reduces to: 



rad : n < —1 



dust : n < 



-iuj + 2) 
2{lu + 1) 



1 



as —> 00 



(48) 
(49) 



where (^) approaches the general relativity value for cj ^ 00, as expected and (48) agrees with the radiation case 
studied above. 
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VII. EXACT SOLUTIONS TO THE FLATNESS PROBLEM 



A. Radiation-dominated case 

The exact solutions found for the VSL theory with c = coa" in ||^ can be generahzed to the BD case if one assumes 
that the variation of the speed of hght is governed by a relation of the form 



c ^ co{ay/G(t)y' (50) 

This reduces to the previously studied (non-BD) case when is constant. In the Einstein frame one then has c = cqo" 
and in the radiation-dominated epoch the first of equations ( p^ ) becomes 



Hence, for radiation, one has the integral 



pa'^^B+ (52) 



AnG{2n + 2 

with B constant. Since pa* = pa*, in the Jordan frame, one has 

This solution can be verified directly from eq. (||), although it would have been difficult to guess it without a foray 
into the Einstein frame. 

Conditions for solving the fiatness problem can now be derived by inspection of the Friedmann-like equation: 

a J <pa'^ (pa b \<p I 

Since (j) approaches an asymptotic value we still require that n < — 1 for the term in Ci to dominate the curvature 
{K) term at late times. 

If (p is decreasing we do not need the speed of light to decrease in time so fast in order to solve the fiatness problem. 
Since 

c a 1(f) 

- = "-+9^ (55) 
c a 2 (p 

we see that while 4>/4> is non-negligible we have \c/c\ < \n\a/a. Weaker gravity in the early universe therefore assists 
a varying speed of light in solving the fiatness problem 



Introduce the variable 



and assume that c varies as 



Hence, (pq) becomes 



Integrating, we have 



B. Other equations of state 

X = 0a3T-2 (56) 
c = coa;"/^ (57) 



OTT 



3i4:c^n(/)"+ia(3'/-2)(n+i) 



if n 7^ — 1. This solution can be used to study the evolution for general 7. 
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VIII. SOLUTIONS TO THE QUASI-FLATNESS PROBLEM 



A. The constant (h case 



We want to discover if it is ever natural to have evolution which asymptotes to a state of expansion with a non- 
critical density (for example, say, flo ~ 0.2, as some observations have implied). For simplicity we consider first the 
solutions with constant (j). The conservation equation with constant (p reduces to 



a^i ~ 47r ' 



2n-3 • 



which integrates to give 



^ flS-/ 47r(37 + 2n-2) ^ ' 



with B constant, so substituting in (Q) we have 



a2 3(/.a3^ a2(i-«)(37 + 2n- 2)' 

iFrom. this we can easily determine the attractors at large a. Specialising to the radiation case, we have 

a? _ 8ttB K cl 

a2 " 30a4 ~ a2(i-n)(^i 

Now, the density parameter is defined by: 

For a quasi- flat open universe the ratio between the two terms on the right-hand side is approximately constant. From 
the solution with 7 = 4/3 we have 



SttG I B I 3Kcln \ 

n _ 3 87rG(n+l)a2(i-'«) J 



(61) 



If n < — 1 the energy production term (second term on the right-hand side of ( |59|)) is subdominant at late times, and 
so this ratio goes to infinite, meaning 17 = 1. IfO<n<— 1 then the second term in (|5^ ) dominates at late times 
and therefore the expansion asymptotes to one displaying 

^ " (62) 



71+ 1 

that is r2 = — n. So the key feature is that in these scenarios the curvature terms associated with violations of energy 
conservation dump energy into the universe at the same rate as the curvature term in Friedmann equation. Therefore 
the ratio of the two terms in the Friedmann equation stays constant, leading to an open universe with finite f2 value 
today. 

B. The varying (j> case 

The formulation of the radiation case in terms of the variable ?/(t) — (f)a? allows us to extend the analysis to the 
varying-(/) case. Taking 

c = coy™ (63) 

we have 
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— AK c y 



and 



47r(2m + 1) 



(64) 



(65) 



We introduce the density parameter 



so that 



n = — 

Pi 



p _ y'^-A^l + ^)+4Kc^y^ 



y''-A^l + ^) 



Stt 



n-1 SKcf^y 



2,,2m+l 



B 



3KmcQy 



,2,,2m+l 



47r(2m + 1) 



Asi— >oo,r^oo,?/^oowe see that f2 ^ 1 if 2m + 1 < (in Section we give a full exact solution for the case 
m = — I which falls into this class), but if 2m + 1 > we have approach to a quasi- flat open universe with 

n -2m 

Again, we can have a natural quasi- flat asymptote with < < 1. 



C. General asymptotic behaviour 



Consider the behaviour of eq. (^9|) at large r and y in the case where 2m -I- 1 > 0; that is, where we have the 
quasi- flat attractor f2 — > fico = —2m as t — > cxd. This assumes the constants in the Friedmann equation allow sufficient 
expansion to occur (so there is no collapse at a finite early time). Since 



we have 



as T ^ 00. Using (EO) we get 



so, since fioo < 1, we have 



Thus t - and 



y -^y 



m+1 



-(r + ro) 



= 00 exp 
t>o as T — > 00. So, we have 
a(T) ~ ■ 



exp 



2(1 -o^) 



a ^ t(i+noo) 



as i — > 00. When Q^o = 1 we have the expected a ^ flat radiation asymptote. 
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IX. THE LAMBDA AND THE QUASI-LAMBDA PROBLEMS 



A. The General Relativity Case 



Let us consider the impact of a varying speed of light on the general relativity case. Similar results will occur 
in the BD case (to which it reduces exactly in the case of radiation with constant 0). If we wish to incorporate a 
cosmological constant term, A, (which we shall assume to be constant) then we can define a vacuum stress obeying 
an equation of state 



where 



PA 



PA 



-PAC 



SttG' 



Then, since G is constant, and replacing p by p + pA in (Q), we have the generalisation 

a p iKcc 
p + 3- p + — +PA = 



We shall assume that the matter obeys an equation of state of the form ( |23| ) and the Friedmann equation is 

a2 _ SttGp Kc^ Ac^ 



(66) 



(67) 



(68) 



(69) 



We also assume that c = coa" again, so eq. ( [Sq ) integrates immediately to give 



p 



47rG(2n - 2 + 37) 47rG(2n + 37) ' 



with B a positive integration constant. Substituting in (69) we have 



_8TrGB Xc^a^f"-!) (2 - 37) Ajcla^" 
— _ I \ 



a 



(37 + 2n) 



(70) 



ia^-t (2n - 2 + 37) 

Eq. (|70| ) allows us to determine what happens at large a. 

If —37 > 2n > 2n — 2 then we see that the flatness and lambda problems are both solved as before. There are three 
distinct cases: 



1. Case 1: 2n > -87 > 2n - 2 



The A term dominates the right-hand side of (70), the curvature term becomes negligible, and 

A7Coa^" 



So, at large t,we have 

Note that for radiation this case requires 

while for dust it requires 

For general fluids it requires 



a2 (37 + 2n) ■ 



-1 < Tlrad < -2, 



2 < ''^dust < 



37 ^ ^ 2-37 

-7:- < n < — - — . 
2 2 



(71) 
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2. Case 2: 2n > 2n - 2 > -87 



The scale factor approaches ( [7l| ) but the curvature term dominates the matter density term. Define 



and then we have that, at large a, 



n 



SttGB 



n-1 



3Kcl 



■a 



3Kc^ 



2-2n-37 



Thus ri/(f2 — 1) ^ 00 as a — > 00 for —87 < — 2n + 2. But when 2n > —87, 



n - 1 



K{2n + S-f) 



00 



in (|70D, and so ^ 1. 
If A = we note that 



2n 



2n - 2 + 87 



when —87 < 2n — 2 and this i s just the solution for the quasi-flatness problem found above for general relativity and 



Brans-Dicke theory in Section VIII 



For A 7^ 0, when the A term dominates at large a we see that 

D,m p 2n 
pa 



2n + 3-f 

and we have a 'solution' to the quasi-lambda problem {ie the problem of why and ^a are of similar order today). 
Recall that in this case we have n < 0, 2n + 87 > and so the asymptote is again of the form a ~ t~ . 



X. CHALLENGES FOR QUASI-FLAT AND QUASI-LAMBDA SCENARIOS 

Scenarios in which the quasi-flatness problem is solved are considerably more exotic than the VSL solution to the 
flatness problem. Unlike flat scenarios they have a Planck epoch, something which may be a problem. We discuss 
this issue in the Appendix. 

These scenarios possess several other unusual features. In standard flatness VSL scenarios, the expansion factor in 
the radiation-dominated phase is still a cx t^^^. Standard nucleosynthesis should still be valid unless there are changes 
to other aspects of relevant strong and weak interaction physics (which seems likely). However in scenarios which 
solve quasi-flatness we have a oc , which for a fl = 0.2 attractor means a (x i"'^^. This could easily conflict with 
the nucleosynthesis constraints. However it is not enough to state that the expansion factor at nucleosynthesis time 
is different: the couplings, masses, decay times, etc, going into nucleosynthesis are all different [|ll|. One must rework 
the whole problem from scratch before ruling out these scenarios on grounds of discordant nucleosynthesis predictions. 

Structure formation is another concern. It was shown in |^ that the comoving density contrast A and gauge- 
invariant velocity v are subject to the equations: 

A' - ^8(7-l)^-h- )A = -7fcw-2— (7-1)Ht (72) 
a 

+ " r- fcc(7-l)(:^ + i£j(T) Inr (78) 

where k is the comoving wave vector of the fluctuations, and F is the entropy production rate, Ht the anisotropic 
stress, and the speed of sound Cs is given by 





—"fkv 




if 




kc{"/ 
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p \ 3j c a' J 



(74) 



Note that the thermodynamical speed of sound is given by — {dp/dp)\s- Since in standard Big Bang models 
evolution is isentropic: = {dp/dp)\s = p/ p = p' I p' ■ When c 7^ the evolution need not be isentropic. However, 
we keep the definition — p' / p' since this is the definition used in perturbative calculations. One must however 
remember that the speed of sound given in ( |7^ ) is not the usual thermodynamical quantity. With this definition one 
has 5p/ 5p = p' / p' for adiabatic perturbations; that is, the ratio between pressure and density fluctuations mimics the 
ratio of its background rate of change. 

Let us assume a radiation-dominated background. For superhorizon modes {ckrj ^ 1) there is a power-law solution, 
A (X rj/^ with f3 = 2(n -|- 1) and f3 = n — 1. The general solution takes the form: 

A = + Btj"-^ (75) 

where A and B are constants in time. For a constant c {n = 0) this reduces to the usual A cx 77^ growing mode, and 
A cx decaying mode. If the flatness problem is to be solved, one must have n < — 1. If this condition is satisfied 
there is no growing mode. 

This is an expression, in the context of Machian BD scenarios, of the link between solving the flatness problem and 
suppressing density fluctuations. Flatness is imposed in VSL by violations of energy conservation, acting so as to leave 
the universe is a state with = 1. This process acts locally, so it also suppresses density fluctuations. Alternatively, 
we can see that the approach to flatness everywhere means that inhomogeneous variations in the spatial curvature 
must also die away. 

Accordingly, we see that in VSL scenarios there is a strong connection between solving the flatness problem, and 
predicting a perfectly homogeneous universe. Quasi-flat scenarios therefore risk not solving the homogeneity problem. 
On the other hand there could some mechanism for amplifying thermal fluctuations to become seeds for the large scale 
structure of the universe jsj. Indeed the equations above should provide a transfer function converting the thermal 
white noise spectrum into a tilted or flat spectrum. If this were the case, then these theories would predict a link 
between the spectral tilt and fl. 

Finally, one may wonder whether such c{t) could ever arise in a dynamical VSL theory. In work in preparation 
psf we address this issue, with the result that indeed if ip = log{c/co) is a scalar field with a Brans Dicke type of 
dynamics: 

^ + 3-Tp ^ AnGujp (76) 
a 

then one has a Machian solution with an exponent n related to the coupling lo of the theory. Hence there will be a 
range of coupling values for which the fiatness problem is solved, and a range for which the quasi-flatness problem is 
solved. In the latter case, not only can we predict an open attractor but also the value of fi of the universe is related 
to a coupling constant. 



XI. CONCLUSION 



We have performed an extensive analysis solutions to Brans-Dicke theories, with varying G, in which the speed 
of light is also permitted to vary in time. We have found cosmological scenarios with novel features. For power- 
law variations in the velocity of light with the cosmological scale factor we identified the cases where the flatness 
problem can be solved. These generalise the conditions found in earlier investigations of this VSL theory. Unlike in 
inflationary universes which solve the flatness problem, no unusual matter fields are required in the early universe. 
We have also identified the cases in which the quasi-fiatness problem can be solved; that is, where there can be 
asymptotic approach at late times to an open universe with a density close to that of the critical value. Similarly, we 
identified those variations of c which provided solutions of the lambda problem and the quasi-lambda problem. The 
possibility of solving the quasi-fiatness and quasi-lambda problems in this way is a genuine novelty of the VSL theory 
that distinguishes from the standard inflationary universe scenario. We have also discussed some problems with the 
VSL scenario, highlighting in particular the matter-radiation transition solutions and the role of the Planck epoch in 
setting initial conditions (see Appendix for further details). We have also examined the behaviour of inhomogeneous 
perturbations to the homogeneous and isotropic solutions and found the conditions for density perturbations to grow 
or decay. 

Our investigations reinforce the conclusions of our earlier investigations of varying-c theories without varying G: 
the scope for obtaining cosmological models which share a number of appealing properties, which closely mirror those 
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of the observed universe, suggest that cosmologies with varying c should be thoroughly explored. It is a challenge 
to find observational predictions which would allow future satellite probes of the microwave background radiation 
structure to distinguish them from inflationary universe models (with constant c) . We hope that this paper will serve 
as a further stimulus to undertake those investigations and to search out new ways of testing the constancy of the 
traditional constants of Nature flfll . 
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APPENDIX - PLANCK TIME IN VSL SCENARIOS 

The Machian VSL scenario, in which c — cqc", introduced by Barrow has significant advantages to the phase 
transition scenario, in which the speed of light changes suddenly from to c^, preferred by Albrecht and Magueijo 

. In the phase transition scenario one runs into the problem of having to decide when to lay down "natural initial 
conditions" (that is e and eA of order 1). For a phase transition occurring at time tc ^ tp, the Planck time tp (built 
from the constants as they were before the transition) is much smaller than tc- But why should we lay down natural 
initial conditions just before the phase transition? If the only scales in the problem are the ones set by the constants 
as they were before the transition, then natural initial conditions should be set at tp. If we are to lay down natural 
initial conditions at tp then the universe goes off the attractor well before the phase transition. A catastrophic phase 
precedes the phase transition, in which the universe becomes Milne (curvature dominated) or dc Sitter (A dominated). 
Albrecht and Magueijo noted that this catastrophic phase is not the end of the universe in VSL scenarios. A varying 
speed of light would not write off Milne or de Sitter universes, but would still push them towards an Einstein-de Sitter 
universe. The only universes which would be selected out in this process are the ones with positive curvature, which 
would end in a Big Crunch, well before the phase transition. 

The Machian scenario does not have this problem, if we are content with solving the flatness but not the quasi- 
flatness problem. If we consider the radiation dominated phase, c — coa" with n < — 1, and a oc then we have 
t/tp oc t^^^. Hence there is no Planck time in these scenarios: as t ^ 0, one has t/tp oo. As we go back in 
time, the universe becomes hotter and hotter, but the Planck temperature also increases, and is never achieved at 
any time. The idea of setting natural initial conditions at Planck time does not make sense in these scenarios. We 
have a universe constantly pushed towards an attractor, which is flat, and has zero cosmological constant. 

Scenarios in which the quasi flatness problem is solved do not have this desirable feature. We find that t/tp cx i^-" . 
Hence we must have a Planck time in our past in these scenarios. 
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